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Abstract 

Let a continuous random process X defined on [0, 1] be (to + /3)-smooth, Q < m, 
< /3 < 1, in quadratic mean for all <: > and have an isolated singularity point at i = 0. 
In addition, let X be locally like a TO-fold integrated /^-fractional Brownian motion for all 
nonsingular points. We consider approximation of X by piecewise Hermite interpolation 
splines with n free knots (i.e., a sampling design^ a mesh). The approximation performance 
is measured by mean errors (e.g., integrated or maximal quadratic mean errors). We 
construct a sequence of sampling designs with asymptotic approximation rate for 
the whole interval. 

Keywords: Approximation, random process, sampling design, Hermite splines 

1 Introduction 

Let a random process X{t), t £ [0, 1], with finite second moment be observed in n points and 
have different quadratic mean (q.m.) smoothness at i = (i.e., an isolated singularity). We 
consider piecewise polynomial approximation of X, combining two different Hermite interpo- 
lation splines for the interval adjacent to the singularity point and for the remaining part (i.e., 
a composite Hermite interpolation spline). A sequence of sampling designs (i.e., meshes) is 
constructed to improve the asymptotic approximation performance (e.g., integrated or maxi- 
mal q.m. errors) eliminating the effect of the singularity point. The proposed technique can 
be applied also to random processes with a finite number of isolated singularity points. In 
principle, the main idea is well-known for nonlinear approximation of deterministic functions 
(see, e.g., de Boor, 1973, DeVore, 1998, and references therein). We use a finer mesh where 
the target function is singular and a coarser mesh where it is smooth. The primary question 
remains however: how to measure this smoothness for random functions in order to obtain 
definitive results. 

More precisely, let l-th derivative of X satisfy an a-Holder condition < Z,0 < a < 1, 
for [0, 1] . It is known that the approximation rate optimal in a certain sense for 

such class of processes and Hermite splines (see, e.g., Buslaev and Seleznjev, 1999, Seleznjev, 
2000). Let, additionally, X have a continuous m-th derivative satisfying a local stationarity 
condition, Herman (1974), with parameter f3, and 0<m,0</3< 1,1 + a < m + 13, for 
all points t G (0,1]. We approximate X by a corresponding composite Hermite interpolation 
spline and set up a sequence of sampling designs attaining the asymptotic approximation rate 
^-{m+p) £qj. ^YiQ whole interval [0, 1]. 
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The nonlinear approximation results obtained in this paper are related to various prob- 
lems in signal processing (e.g., optimization of compressing digitized signals, see, e.g., Cohen 
and D'Ales, 1997, Cohen et al., 2002), in numerical analysis of random functions (see, e.g., 
Benhenni and Cambanis, 1992, Kon and Plaskota, 2005, Creutzig and Lifshits, 2006, Creutzig 
et al., 2007), in simulation studies with controlled accuracy for functionals on realizations of 
random processes (see, e.g., Eplett, 1986, Abramowicz and Seleznjev, 2008). Hermite spline in- 
terpolation for continuous and smooth random functions is studied in Seleznjev (1996, 2000), 
and Huesler et al. (2003). Ritter (2000) contains a very detailed survey of various random 
function approximation problems. 

The paper is organized as follows. First we introduce a basic notation. In Section[2l we con- 
sider composite spline interpolation for certain Holder type classes of random processes, which 
behave locally like m-fold integrated fractional Brownian motion for all intervals [a, b] C (0, 1] 
and construct sequences of sampling designs with asymptotically optimal approximation rate. 
In the second part of the section, the approximation of more general Holder's classes of random 
functions and the accuracy for the corresponding composite Hermite spline approximation are 
studied. Sections [3] and H] contain the results of numerical experiments and the proofs of the 
statements from Section O respectively. 

1.1 Basic notation 

Let X = X{t), t G [0, 1], be defined on a probability space {VL,^,P). Assume that for every 
t, the random variable X{t) lies in the normed linear space L^(r2) = Lp'{Q.,^ ,P) of random 
variables with finite second moments and identified equivalent elements with respect to P. We 
set II Z II := (EZ2)V2 for all Z G L'^iVL). Let C'^dO, 1]) denote the space of random processes 
with continuous q.m. derivatives up to order m > 0, and C""([0, 1]) be the corresponding space 
of non-random functions. Let "PfedO, 1]) be the space of stochastic polynomials of order A; > 0. 
We define the norm for any X £ ^"([0, 1]) by setting 



oo, then 1/p = 0. Denote by f{t) x g{t) as i — )• a the weak equivalence property, i.e., 
Ci git) < fit) < C2 g{t) for some positive Ci, 6*2,6 and for all t G {a,b] and similarly for 
On ^ 6„ as n — )■ 00. 

Let X be sampled at the distinct design points T„ := (to, ii, • • • , in) (also referred to as 
knots), and the set of all (n + l)-point designs be denoted by Dn ■= {Tn : = to < ti < ■ ■ ■ < 
t„ = 1}. We suppress the argument n for the design points tk = tkin) from T„, k = 0, . . . ,n, 
when doing so causes no confusion. Recall that, for any / G C^dO, 1]), m > 0, the piecewise 
Hermite polynomial Hk{t) = ff/j(/, T„)(t), of degree fc = 2r + 1, < r < m, is the unique 
solution of the interpolation problem H^\ti) = f^^\ti), where i = 0, ...,n, j = 0, ...,r. 
Analogously, we suppose that for X G C"^([0, 1]), the process and its first r < m derivatives 




1 < p < 00 



and 1 1 X 



00 



max[o,i] II Xit) II for p = 00. Henceforth, we use the convention, that if p 
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can be sampled, and write Hk{X,Tn) with k = 2r + 1 < 2m, + 1 to denote a corresponding 
stochastic Hermite sphne. Define Hq^j^{X,Tn), g < /e, to be a composite Hermite spline 



Hq{X,Tn)it), tG[0,ti] 

Hk{x,Tn){t), te[ti,i] ■ 



To formulate the asymptotic results, we introduce a quasi regular sequence (qRS) of sampling 
designs {T„ = T„(/i)} generated by a positive density function h{t), t G (0, 1], via 



/ h{t)dt = i/n, i = l,...,n, / h{t)dt = 1. 
Jo Jo 



(1) 



Henceforth we assume that h{t) is continuous for t G (0, 1] and if h{t) is unbounded in t = 0, 
then h{t) — >■ +oo as t — >■ 0+, otherwise h G C([0, 1]). We denote this property of 
by: {Tn} is qRS{h). If the function h is positive and continuous on the whole interval [0, 1], 
the corresponding generated sampling designs are called regular sequences, RS(/i), and were 
introduced by Sacks and Ylvisaker (1966) for certain time series models. In particular, if h 
is uniform over [0, 1] {h{t) = 1, t G [0, 1]), then the regular sampling becomes the equidistant 
sampling including the cndpoints. For random process linear approximation problems, see Su 
and Cambanis (1993), Seleznjev (2000), and references therein. Define the related distribution 
functions ^ ^ 

H{t):= [ h{u)du, G{t) := H'^it) = [ g{v)dv, te[0,l], 
Jo Jo 

i.e., G is a quantile function for the distribution H. Then by the definition, 

ti = G{i/n), i = l,...,n, 

the knots are i/n-quantile points of H. The quantile density function 

g{s) := G'{s) = l/h{G{s)) 

is assumed to be continuous for s G [0, 1] with the convention that g{Q) = if h{t) — )• +oo as 
t — )■ +0. Following Sacks and Ylvisaker (1966), we define asymptotic optimality of a sequence 
of sampling designs T* by 

lim ||X - Hq,k{X,T:)\\J mi \\X - Hq,k{X,T)\\p = 1. 

Now we introduce the classes of processes used throughout the paper. We say that: 

(i) X G C™''^([a, b], V{-)) if X G C™([a, b]) and is locally Holder continuous, i.e., if for all 

t,t + s,e [a,b], 

||X("^)(t + s)_X(™)(t)|| <V{t)'^/^\sf,0 < p <1, (2) 

for a positive continuous function V{t),t G [a,b], and some t G [t,t + s]. In particular, if 
V{t) = C, t G [O'jb], where C is a positive constant, then X^"^^ is Holder continuous, and we 
denote it by X G C"^'^([a, b],C) 



3 



(ii) X e B'^'^{[a, b],c{-)) if X G C"^{[a, b]) and there exist < /3 < 1 and a positive continuous 
function c{t), t G [o.,b], such that X^"^^ is locally stationary, i.e., 

1 1 v{m) (f I a\ _ vim) (f \ 1 1 

lim ^ \s\^ " ''^^^^^^ uniformly in t £ [a, b], (3) 

(iii) X G C;S'^'''((0, 1], c(-), if there exist < /3 < 1 and positive continuous functions 
c(t), t G (0, 1], such that X G C™'^([a, 6], n S™'^([a, 6], c(-)) for any [a, b] C (0, 1]. 

By definition, we have that V{t) > c{t),t G (0,1]. Moreover, if X G C("'+i)((0, 1]), then 
c(t) = II X^"^'^^\t) Ip and we may set also V{t) = c{t), t G (0, 1]. Some properties of locally 
stationary processes are considered in Berman (1974), Hiisler (1995), and Seleznjev (2000). 

Example 1 Let B{t),t G [0, 1], be a fractional Brownian motion with the covariance function 
r{t,s) = (|tp'^ + |sp'^ ~ |i ~ where < /3 < 1 is a Hurst parameter. Define a time 

changed version of the process B, 

X{t) ■.= B{Vi),te [0,1]. 

Then 

\\X{t + s)-X{t)\ \ = \VtTs-Vif < |s|^/2 for alH,t + s G [0,1]. 
However on the nonsingular part of the domain, we get 

\\X{t + s)-X{t)\\ = {VtTs + Vi)-^\sf for alH,t + s G (0,1] 

and therefore X G C^'f^/'^{[0, 1], 1) n CB^'^i{0, 1], c(-), V{-)) with V{t) = c{t) = (4t)-^. 

For the differentiable case, a representation of the approximation error can be obtained as 
a consequence of the following proposition (see Seleznjev, 2000), which is a q.m. variant of the 
well-known Peano kernel theorem (see, e.g., Davis, 1975, p. 69). We formulate this proposition 
for further references for a closed interval [a, b] C [0, 1]. Define linear operators of the following 
type over C'^{[a,b]), 



(4) 



where the integrals and derivatives are taken in quadratic mean. The functions dj(-) are 
assumed to be piecewise continuous over [a, 6] and all points Sij G [a, 6]. Further, for any 
s G [a, 6] and a given m, let (i)+ := max(0,t), and Ps,m{t) := (t — s)™, t G [a,b], be the 
function of t for a given s. 

Proposition 1 (q.m. Peano kernel theorem) Let R be of the form ^ and R(Y) = if 
Y GVrr,{[a,b]). Then for all X G C""+H[o,6]), 

RiX)= t X^-^+^\s)K„Xs) ds, (5) 



where the Peano kernel Km{s) := R{ps,m) /fnl, s G [a,b]. 



4 



Example 2 For a given point t in [0, 1], denote by Km,kit, s) the Peano kernel (as a function 
of s) of Hermite interpolation of X G C"^~^^([a,b]) by H^it). In particular, for the two-point 
piecewise linear interpolation of X G C^([0,1]) on the interval [0,1], we have KQ^i{t,s) = 
{t - s)+-t, t,se [0,1]. 

For a zero mean fractional Brownian motion B{t), t G [0,1], with Hurst parameter /3 the m-fold 
integrated Brownian motion Emit) is defined by setting 

1 

Bm{t) := 7 TTT / {t - s)'^-'B{s) ds, m>l, Bo{t) = B{t), (6) 

(m-l)!7o 

S^)(0) = 0, J = 0, ...,m - 1, BiT\t) = B{t). Let Rk,2{Bm){t) be the remainder for the 
two-point Hermite interpolation of B„i{t), t G [0, 1], by the Hermite spline Hk{Bm){t) with the 
norm 

•= \ \Rk,2{Bm) Hp- 
Explicit expressions for particular values of h^'^ can be found in Seleznjev (2000). 

Recall that a positive function / is called regularly varying (on the right) at the origin with 
index p if for all A > 0, 

4^ ^ A'' as X ^ + . (7) 
/(^) 

We denote this property by / G lZp{{)+). Throughout the proofs we use the following properties 

of regularly varying functions (see, e.g., Bingham et al., 1987): 

(Rl) convergence in ([7]) is uniform for all intervals 0<a<A<6<oo; 

(R2) if f{x) G 7^p(0+) and F{x) = f{v)dv, x G [0, 1], then F G 7lp+i(0+). 

Write / G 7l+(r(-), 0+), if for some r G 7lp(0+), a > 0, we have f{t) < r{t) for all t G (0,a]. 

2 Results 

We consider first the case when an optimal approximation rate can be attained by using a 
suitable composite splines and sequences of sampling designs (or generating densities). Let 
X G C''"([0,1],M) n C^"^'^((0,l],c(-),V^(-)), I + a < m + p. We formulate the following 
condition for a local Holder function V and a sequence generating density h: 
(C) let g G 7^+(r(•),0+), where 

r{s) = o(s('"+/3)/('+°+i/p)-i) as s ^ 0; (8) 
if p = oo, then y(t)i/2r(i?(t))'"+^ ^ as t ^ 0; 

if 1 < p < oo and, additionally, V{G{-)Y/'^ G 7^+(i^(•), 0+), then R{H{t))r{H{t))"'+l^ G Lp[0,b] 
for some b > 0. 

In the following theorem, we describe the class of generating densities eliminating the effect of 
the singularity point for the asymptotic approximation accuracy. 
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Theorem 1 Let X £ C''"([0, 1], M) r\CB'^'l^{{0,l],c{-),V{-)), l + a <m + P, with the mean 
f G C™'^([0, 1], C), f3 < 6 < 1, be interpolated by a composite Hermite spline Hq j^{X,Tn), 
I < q, < k, where Tn is a qRS(h). Let for the density h and the local Holder function V , 
the condition (C) hold. Then 

Jim n^+'^llX - i7,,fc(X,r„,)||p = h'^;^^\\c"^h-^^^P\\^ > 0. (9) 

Remark 1 (i) The condition (C) is technical and we conjecture that the main resuh of the 
above theorem is vahd for a wider class of processes and generating functions. In particular, 
for p = oo, if in (C) the following condition is used, 

V{t)^'^r{H{t)r+^ ^Dast^O (10) 

for some constant Z? > 0, it follows straightforwardly from the proof that 

||X-/7g,fc(X,r„)||oo ><n-("+^) as n^oo. (11) 

More general conditions are considered also in Proposition [2j 
(ii) Let 7 := l/(m + /3 + 1/p) and write 



h*{t) := c(t)^/y c{sy/^ds, t G (0, 1]. 



It is well known that h* minimizes the asymptotic constant in ([9]) (see, e.g., Seleznjev, 2000). 
Thus if the conditions of the above theorem hold for X and h* , then h* is the asymptotically 
optimal density and 



lim n^+^WX - H,^k{XMh*mp = 6™;f \\c'/'\ 



7- 



For example, if c{t) ~ Ct for some constant C > and V{t) x c{t) as t ^ 0, then for any 
1 < p < oo the inequality 

m + p J I 7 

is a sufficient condition for h* to be the asymptotically optimal density, as it implies the condi- 
tion (C). However for p = oo, only (llOp holds and it follows from (jlip that the approximation 
error is of order for large n. 

Example 3 Consider X{t),t G [0, 1], defined in Example 1. By previous calculations we have 
c{t) = (4t)-^,0 < /3 < 1, and 7 = l/{(3 + 1/p). Now the condition ([H]), i.e., 

is satisfied for any 1 < p < oo, hence 

r(t) = (1-/37/2)^-^^/2, tG(0,l], 
is the asymptotically optimal density. 



6 



We proceed to some cases not included in Theorem [T] but important for applications. Specif- 
ically, like in the conventional approximation for interpolation by piecewise polynomials (or 
splines) of order k and smoothness m > k, the upper bound for the optimal approximation rate 
is n~^''^^\ We investigate the corresponding generating densities leading to the asymptotically 
optimal solution in such setting. Namely, let X G C''"([0, 1], M) n ^"((0, 1]) be approximated 
by a composite Hermite spline Hg^^- For any q < k < m we have Cfc(i) := ||X*^'^"*'^''(t)|p. We 
introduce the following modification of the condition (C), where /c + l is used instead of m-|-/3: 
(C) let g G 7^+(r(•),0+), where 

r{s) = as s ^ 0; (13) 

if p = oo, then Ck{t)^/'^r{H{t))^+^ ^ as t ^ 0; 

if 1 < p < oo and, additionally, Cfc(G(-))^/^ G 7^+(i^(•), 0+), then R{H{t))r{H{t))''+^ G Lp[0,b] 
for some 6 > 0. 

Theorem 2 LetX e C''"([0, 1], Af)nC"^((0, 1]), l+a < m, with the mean f e C7"^+i'^([0, 1], C), 
<6 <1, be interpolated by a composite Hermite spline Hq^f^{X,Tn), I < q, k < m, where Tn 
is a qRS(h). Let for the density h and the function V, the condition (C') hold. Then 

lim n^+i||X - F,,fc(X,r„)||p = bl^l WcJ'h-^'^+'X > 0, (14) 

where b^^l = B{p{k+l) /2+l,p{k+l) / {k+l)\, l<p<oo,and bl'^^ = 2-(^+i)/(A;+l)!, 
and B{-,-) denotes the beta function. 

Similarly to the above case m < k, we consider an optimal density for such spline approxima- 
tion. 

Remark 2 Let 7^ := l/{k + 1 + 1/p) and write 

h*{t) := Cfc(t)^^-/y Ckisy^/^ds, t G (0, 1]. 

If the conditions of Theorem [2] hold for X and h* , then h* is an asymptotically optimal density 
and 

lim n^+'\\X - H,,k{XMh*mp = blf^ ||cf ||^,. 



The following proposition says that an improvement of the approximation rate is gradual 
with respect to some properties of knot generating densities. We introduce conditions for 
a more general class of processes and generating densities when approximation rates n~'^, 
K £ [I + a,m + (3], are attained. 

Proposition 2 Let X e C''"([0, 1], Af) n C"'''^((0, 1],V{-)), l + a <m + f3, with nonincreasing 
V{-), be interpolated by a composite Hermite spline Hq^k{X,Tn), I < q, m < k, for sampling 
designs Tn such that 

hi < Af~^/('+''+^/P)n~''/('+"+^/P\ 

h, < y(t^._,)-i/(2(™+/3+iM)^-'^/(-+/3+i/p)^ i = 2, . . . , n, ^^^^ 
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where k £ [I + a,m + f3]. Then 



n^\\X -Hg^kiX,Tn)\\p<C 

for some positive C . 

Remark 3 The constant C in the above proposition depends only on the corresponding Peano 
kernels of the composite spline and hence some conventional (deterministic) approximation 
theory results can be used (see, e.g., Davis, 1975, p. 69). 

3 Numerical experiments 

In this section, we present some examples illustrating the obtained results. For given knot 
densities and covariance functions, first the pointwise approximation errors are found analyt- 
ically. Then numerical maximization and integration are used to evaluate the approximation 
errors for the whole interval. Notice that the nonlinear approximation here gives a significant 
gain in the order of approximation error. We consider the following knot densities 

hx{t) = ^t-^-\ tG(0,l], A>0, 

say, power densities. This leads to the sampling points ti = (i/n)^ ,i = 0, . . . , n, where A = 1 
corresponds to the uniform knot distribution (equidistant sampling). Define by 

6n{X,Hg^k,hm ■.= X{t)-H^^k{X,Tn{h)){t), t G [0,1], 

the deviation process for the approximation of X by the composite spline Hq^k with n knots 
generated by the density h and write 

e„(/i) = en{h,X,Hq^k,p) ■= \ \Sn{X, H^^k, h)\\p 
for the corresponding mean error with the density function h. 

Example 4 Consider X{t),t £ [0,1], defined in Example 1 and let (3 = 0.4. Let en{hx) = 
en{hx, X, Hi^i,oo) be the mean maximal approximation error, p = oo, when the both splines 
are piecewise linear interpolators, q = k = 1. For the power density h\ and X, the conditions 
of Theorem [1] are satisfied if A > 2. We choose the following values of the power density 
parameter: Ai = 1 (uniform knot distribution) and A2 = 2.1. Figure a) shows the (fitted) 
plots and evaluated values of the mean maximal errors en{hxi),i = 1,2, versus n (in a log- 
log scale). These plots correspond to the following asymptotic behavior of the approximation 
errors: 

enihx,) ~ Ci ~ 0.377, 

e-n{hx2) ~ C2 n"°-s, C2 ~ 0.295 as n 00. 

For example, the minimal number of observations needed to obtain the accuracy 0.01 is ap- 
proximately 8727 for the equidistant sampling density hx^, whereas it needs only 69 knots 
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Figure 1: (a) The (fitted) plots of e„(/iAi) (daslied-line) and en{h\^) (solid line) versus n in 
a log-log scale, (b) The convergence of nP'^enihx^) (solid line) to the asymptotic constant 
(dashed line). 

when hx^ is used, i.e., Theorem [1] is applicable. Figure [TJb) demonstrates the convergence of 
the scaled approximation error n ' c^ii^hx^^ to the asymptotic constant obtained in Theorem [Tl 

Example 5 Let Y{t), t G [0, 1], be a zero mean stationary Gaussian process with the covariance 
function r{s,t) = exp{— (s — t)^}. We consider a distorted version of the process, 



Then X £ C0'°-9([0,l],l) n C"((0, 1], c™(-)) for any m > 0, where c„(t) = ||X(™+i)(t)||2, 



t G (0, 1]. The process Y has infinitely many q.m. derivatives, hence the approximation rate 
of Hermite spline approximation is limited by the order of spline only. However, the linear 
methods applied to X would suffer a substantial loss of efficiency. Consider now an approx- 
imation of X by the composite Hermite spline Hi 3, i.e., a linear function on the interval 
adjacent to singularity, and the cubic Hermite spline otherwise. We investigate the mean in- 
tegrated {p = 2) approximation error en{h\) = en{hx, X, Hi^3,2). For the corresponding local 
stationarity function, we have 



where C ~ 0.0432. For the power density hx and X, the conditions of Theorem [2] are satisfied 
if A > 20/7. We choose A2 = 3, A3 = 4, and A4 = 5 satisfying this condition, together with 
Ai = 1, the equidistant sampling density. Figure [2]J^ a) shows the (fitted) plots and evaluated 
values of mean integrated errors e„(/iAj, i = 1,2,3,4, versus n (in a log-log scale). 
Note that the estimated approximation rate for the uniform knot distribution is n~^'^, whereas 
the power densities with parameters A2, A3, and A4 result in the regular (nonsingular) case rate, 
namely, n~^. By Remark 2, the asymptotically optimal density 



Moreover, h*{t) ~ hx{t) with A = 45/14 ~ 3.21 as t — t- 0. Figure [2]^b) demonstrates the ratios 
&nihxi)/enih*), i = 2,3,4 and an essential benefit in the asymptotic constant for the optimal 
density h*. 



X{t) ■.= f-'^Y{t),t£ [0,1]. 



C3(t) ~ Ct 



.-6.2 



as t — ^ 
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Figure 2: (a) The (fitted) plots of en{hx-), i = 1,2,3,4 versus n in a log-log scale, (b) The 
(fitted) plots of ratios en{hx-)/en{h*), i = 2,3,4, for the asymptotically optimal density h*. 

4 Proofs 

Proof of TheoremUl We investigate the asymptotic behavior of the q.m. error e„(t) := 1 1 X{t) — 
Hq^k{X, Tn){t) II for any t G tj], i = 1, . . . , n, when the number of knots n tends to infinity. 

An explicit formula when / = or m = is used, and the q.m. Peano kernel representation 
is applied otherwise. We show that the densities satisfying the conditions of the theorem 
lead asymptotically to eliminating the effect of the singularity point. Further, we find the 
asymptotic form of || X — ffg^fc(X, T„) ||p for any density h{-) satisfying the conditions. 

Throughout the proof we utilize the next property of the Peano kernel (see, e.g., Seleznjev 
2000): 

Let Km-i,k{t, s),t, s £ [0, 1], be a Peano kernel for the two-point Hermite interpolation, m > 1. 
li m < k < 2m + 1, then 

/ K:,,,^,{t,s)ds = 0. (16) 
Jo 

Additionally, the following (shifting) property of a regular varying function r is used: 
there exist positive C, /3 such that 

1 s 

r(s) < Cr(v) for all - < - < 2, t; / 0, max(s, v) < /3, (17) 

2 V 

which follows directly from the regular variation property (Rl). 
Note that ([8]) and the regular variation property (R2) imply 

G{s) = a as s ^ 0. (18) 

Now we come to the detailed presentation of the proof. We investigate the asymptotic 
behaviour of the q.m. error e„(t) := || X{t) — Hq^k{X,Tn)(t) \\ for any t G i = 1, . . . ,n, 

when the number of knots n tends to infinity. Consider first the case p = oo, the mean maximal 
norm. Let for a fixed 6 > 0, Js denote the largest integer such that tjg < 5. Then the error 
can be decomposed as follows, 

\\X - Hq^kiX,Tn)\\oo = maxen(t) = max(e„j, j = 1, . . . ,n) = max(S'i, ^2, ^s). 
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where e„j := max[t^_^^tj] en{t) 



5i 



52 

S3 



Si{n) = e„,i, 



S2{n) = max(e„j, J = 2,..., Js), 

83(71) = max(e„j,j = J5 + 1, . . . , n). 

For the first mterval and Si, if / = 0, then a piecewise linear interpolator Hi{X,Tn) is 
apphed. Hence, by the definition, (R2), (C), and (fTSj) . we obtain that 



= n'^+^o(n-(™+^)) = 0(1) as n ^ 00. 



(19) 



If / > 0, Proposition [U and (dH) imply, that for t G [0,ii], 

X{t) - Hg{X,Tnm = h{ l\z{s) - Z{0))KUJi,s) ds, 

Jo 

where Z{s) = + hjs), and -f^/L^ ^{t, s) denotes the corresponding Peano kernel for the 

two-point Hermite interpolation at [0, 1]. Now the Holder property and (fT8|) give 

n'^+^ei^n < 071""+^ h[+'' = Cn'^+^G(l/n)'+" = o (1) as n ^ 00. (20) 

For S2, let first m > 1. Then Proposition [1] implies 

X{t) - Hg^k{X,Tnm = X{t) - Hk{X,Tnm = r X(™)(s)i^„_i,fc(t,s)ds 

= hf CY{s)K*m-i,k(trs) ds 
Jo 



(21) 



where Y{s) = X^"^\tj-i + hjs), and K^_^ ^(t, s) denotes the corresponding Peano kernel for 
the two-point Hermite interpolation at [0, 1]. Using (jl6p . we get 

X{t) - Hg,k{X,Tnm = hf f {Y{s) - y(0)) K*^_^^^{i,s) ds. (22) 

Jo 

Now the Cauchy- Schwartz inequality together with ([2|) give 

en,, < ChJ+^V{wjf'\ (23) 
where V{wjY/'^ = max[t^_^_i^] V{tY/'^,Wj = G{sj), sj G [(j - l)/n, j/n], and 

C = C{Km-i,k) = max | |i^m-i,fc(i, Ol b- 

te[o,i] 

Notice that by definition and the integral mean value theorem for a positive Ci, 

hj = g(^-\-G f = g{vj)- < Ci-, i = 1, . . . , n, (24) 

Vn/ \ n J n n 
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and, therefore, from (j22]) we obtain, Vj £ [{j — l)/n,j/n], 

< CV{wj)^/^r{vjr+l^. 
Now the (shifting) property pT|) imphes for some positive C2, 

n"^+^e„j < CV{G{sj))'/\{v,r+^ 

< C2V{G{s,))'/^r{s,r+^, j = 2,...,Js, (25) 

for small enough 6 by the definition, since Sj,Vj G [{j — l)/n,j/n] and 

1 7 — 1 S-i ?' 

- < < — < — — < 2, j > 2. 

2 j Vj j -I 

For m = 0, we have k = 1, and the explicit expression for the piecewise linear interpolator 
yields 

e„,,(t) = \\X{t) - {l-t)X{t,^,)-tX{t,)\\ < \\Xit)-X{t,^^)\\+t\\X{tJ)-X{t, - 1)|| 



and therefore, by ()24p . the Holder condition, and the shifting property (I17p we get 

n^e„,, < 25(t;,)V(G(s,))V2 < Cr(5,)^l^(G(s,))^/^ (26) 

for some positive constant C. Thus, by (C), (j25p . and (j26p for any e > 0, sufficiently small 
and n large enough, we get 



n 



^+(^82 = max(n-+'^e„,„ J = 2, . . . , Js) < e, (27) 



for any m > 0. 

For S'a, we show first that 

n"^+^53 ^ max(c(t)i/2;^-(™+^) (t)) as n ^ 00. (28) 

The main steps of the proof repeat those of the corresponding result for a regular sequence 
of designs in Seleznjev (2000) for interpolation of a smooth random function by the Hermite 
spline Hk{X, T„). Applying the local stationarity condition and uniform continuity of positive 
c{t),h{t),t G [5,1], we obtain 

n-+/^53 =n™+'' max(e„(t), t G [tj^i,tj],j = J5 + 1, . . . , n) 

=6^f max(c(u;,)i/2/i(u;,)-("+/^),j = Js + l,.. . ,n){l + o(l)), ^^^^ 

for some Wj G [tj-i,tj],j = + 1, . . . , n, and consequently 

n^'+^Ss =6^;f max(c(t)i/2/i(t)-(-+/^),t G [tj,, 1])(1 + o(l)) as n ^ 00. (30) 
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Moreover, 

ma^{c{t)'/'h{t)-(^+^\te[tj„l]) 

= max(c(t)i/2/^(t)-(™+'^),t e [5,1])(1 +o(l)) as oo, 
since yields \tjg — 6\ < hjg < C/n for a positive constant C. 

Let as := max{c{t)^/^h{t)-^"^+l^\t G [tj^, 1]). Note that it follows from (C) that as is bounded, 
and the monotone convergence gives 



as 



t a := sup(c(t)i/2;^(t)-('"+^) , t G [0, 1]) as (5 ^ 0. (32) 



So, first, we select 6 sufficiently small and apply (p7|) . ([30]) . and ([32|) for sufficiently large n. 
Then for the selected (5, (fT9]) and (pOj) imply the assertion. This completes the proof for p = oo. 



The proof for 1 < p < oo is analogous to the previous case, so we give the main steps only. 
Let for a fixed 6 > 0, 

el := / enit fdt = / e„(t)^dt = e^, + <j + E <3 = Si + S2 + S^, 

where 5i = Si{n) := 6^,1, the sum ^3 = S^{n) includes all terms e^j such that [tj_i,tj] C [5, 1], 
say, j > Js + 1, and S2 = - Si - S3. 

For 5i, the first interval, if / = 0, Hi{X,Tn) is applied, hence by the definition, (R2), (C), 
and (jlSp . we have 

„p("^+/3)(e^^^(i))P <^p(m+i)j\^2-P° 2/i?"(it < M2i-P"nP('"+'^)G(l/n)P"+i 

Jo (33) 
=^p("i+/3)Q(^-p(m+/3)) = 0(1) as n ^ 00. 

If Z > 0, Proposition [U and ([16]) imply that for t G [0,ti], 

X{t) - H,{X,Tn){t) = h[ f\z{s) - Z{0))Kl_,^^it,s) ds 

Jo 

where Z{s) = X^^\tj^i + hjs), and i^/L^ ^(t, s) denotes the corresponding Peano kernel for the 
two-point Hermite interpolation at [0, 1]. Now the Holder inequality and (jlSp give 

nP(™+^)(e„,i)P < CnP(™+^)/i}+P'+P" = CnP('"+^)G(l/n)i+P'+P" = o(l). (34) 

For S2 and m > 0, i.e., j = 2, . . . , Js, ([22]) and ([M]) yield that for some positive constant C 

n(™+^)e„,, < Cn(™+^)y(u;,)'/'/i7+^+'/^ < CF(t^;,-)'/'5(^,)^'"+^^/if ^ 

where Wj,Vj £ [{j — l)/n, j /n]. The (shifting) property ([17p together with condition (C) imply 
that for some positive Ci, 

< Ci R{H{t))Pr{H{t))^"'+''^Pdt, j = 2,...,Js, (35) 
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where for Sj = H{uj) G [{j — l)/n,j/n], 

RisMsjr^^ = R{H{uj))r{H{uj)r+^ = mm R{H{t))r{H{t)r+^. 

Thus, for any ei > 0, sufficiently small 6, and large enough n, we obtain 

nim+l3)PsP < Ci f '' R{H{t))Pr{H{t)Y"'+'^^Pdt < ei (36) 

Analogously to the case p = co, the explicit form of piecewise linear interpolator is used to 
obtain (136p when m = 0. 
For S3, we show that 

^(m+/3)p^P ^ (Cp'^)^ / c(t)P/2/j(t)-(™+/3)P(it as n ^ 00. (37) 



Applying the local stationarity condition and uniform continuity of positive c(t), h{t),t E [6, 1], 
we get 

for some Wj E tj], j = + 1, . . . , n, and consequently 

1 



^(m+/3)p^p ^ (^^mj3^f \^ c(t)P/^h{tr^"'+''^Pdtj (1 + o(l)) 

= (C)' c(t)^/2Mt)-(™+/^)^'dt) (1 + o(l)) 

=: 05(1 + 0(1)) as n — )• oo. (38) 

Now it follows by (C), that c(t)P/2/i(t)"(™+^)P is integrable on [0, 1], hence 

asta:= (b'^fY J c{tyP/^h{t)-'^"'+^^Pdt, t E [0, 1]) as <5 ^ 0. (39) 
This completes the proof. 

Proof of Theorem O The proof repeats the steps of that of Theorem [1] for the differentiable 
case. For the calculations of the asymptotic constants we refer to Seleznjev (2000). 

Proof of Proposition 0. The arguments are similar for the first interval and the other intervals 
in both the mean integrated, and maximal norms. So we demonstrate the main steps for the 
[tj^i,tj],j > 2, the mean maximal norm, and the differentiable case. Applying Proposition [1] 
and (fTUj) to the Hermite spline approximation, we have, t E [tj-i,tj],j > 2, 

X{t) - Hk{X){t) = hf [' {Y{s) - Y{0)) K*^.kit, s) ds, 
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where Y{s) = + hjs), and K^__^ denotes the corresponding Peano kernel for 

the two-point Hermite interpolation at [0,1]. Now the Cauchy- Schwartz inequality, and 
the monotonicity of V give 

\\X{t)-Hk{X){t)\\^ < ChJ^^V{t,.if'\ 

where C = C(-R'm-i,fc) = maX(g[o,i] ll-f^m-i,fc(i) Olb- Finally, applying the condition ([15]) im- 
plies the result. 
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